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1. Introduction 

The study of neutral pion production in proton-proton collision near threshold, pp — > 
ppn^, has been inspired by precise measurements of the ncar-thicshold cross section [1, 2]. 
Surprisingly, the measured cross section turned out to be ~ 5 times larger than the early 
theoretical predictions [3, 4]. Main theoretical difficulties in describing the threshold cross 
section may stem from the suppression mechanisms in both one- and two-body 5'-wave 
neutral pion production operators and the typical large momentum transfer k between 
two protons at the threshold, k ~ ^Jm-^nij'^ {k~^ ~ 0.55 fm), where ttItt and mjv are 
the pion and nucleon mass, respectively. Subsequently, some mechanisms to account for 
the threshold experimental data have been suggested in model calculations: one is the 
short-range effect of heavy-meson exchanges [5] and another is the off-shell effect of TriV 
(S-wave isoscalar amplitude in the one-pion exchange production operator [6]. 

Heavy-baryon chiral perturbation theory (HBxPT) is a low-energy effective field theory 
(EFT) of QCD and provides us a systematic perturbation scheme in terms of Q/A^ where 
Q denotes small external momentum and/or symmetry breaking term itIt^ and denotes 
the chiral scale A^ = 47r/7r — 1 GeV: /tt is the pion decay constant. (For reviews, see, 
e.g., Refs. [7, 8, 9, 10, 11].) In Refs. [12, 13, 14] the total cross section of pp — > pp7r° 
was calculated at the tree-level, where the production operators were derived from the 
HB^PT Lagrangian up to next-to leading order (NLO) with the Weinberg's counting 
rules [15], and the matrix element was obtained by DWBA: thus the operators were 
sandwiched between the initial and final two-nucleon wavefunctions obtained by solving 
the Schrodinger equation with the accurate phenomenological NN potentials. Up to one- 
loop order, i.e., next-to- next-to leading order (NNLO), the production operators were 
derived in Ref. [16] and the total cross section was obtained by including the initial and 
final state interactions in Refs. [17, 18]. The other approaches, e.g., a tree-level calculation 
including heavy-mesons [19] and a relativistic calculation including a part of the one-loop 
diagrams [20] were also reported. For a recent review, see Ref. [21] and references therein. 

Though many works on the pp — > pp7r° reaction near threshold in HB^PT have been 
done, some issues in theoretically describing the process have not been fully clarified. 
In the NLO HB^PT calculations [12, 13, 14], a significant enhancement of the off-shell 
TTTrNN vertex function obtained from the NLO HBxPT Lagrangian is found. However, the 
two-body (one-pion-exchange) matrix element with the off-shell ttttNN vertex is almost 
exactly canceled with the one-body matrix element. Thus the experimental data cannot 
be reproduced in the NLO calculations. In the NNLO HB^PT calculations [16, 17, 18], 
a significant contribution comes out of the NNLO corrections and a moderate agreement 
with the experimental data is obtained [17]. However, the chiral series based on the stan- 
dard Weinberg's counting rules shows poor convergence. A modification of the original 
Weinberg's counting rules to account for the large momentum transfer, k ~ ^/Wi^i^rnN^ 
is discussed in Ref. [13]. The production operators at NLO using the modified count- 
ing rules are estimated, and it was reported that the NLO contributions exactly cancel 
among themselves [22]. Recently, some detailed issues for the loop calculations, such as a 
concept of reducibility [23] , a representation invariance of the chiral fields among the loop 
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diagrams, and a proper choice of the heavy- nucleon propagator [24], were also studied.^ 

The problem in those NNLO HB^PT calculations for the process is that no hierarchy in 
constructing the production operators is found. A plausible explanation for this situation 
is that along with the suppression of the leading order Weinberg- Tomozawa term in the 
TTTiNN vertex for the one-pion exchange contribution, the one pion exchange propagator, 
which is usually counted as Q^^, gains another suppression factor Q/niN because of the 
typical large momentum transfer fc^ ~ uiT^mN and thus is counted as ~ {'mT^m]sf)~^ ■ 
Therefore the one-pion-cxchangc diagram for the process could gain the two suppression 
factors (Q/A^)^ and thus become the same order as the one-loop diagrams in the original 
Weinberg's counting rules. Apparently the contributions of the intermediate two-pion 
exchange and the short-range contact diagrams play almost the same role as that of the 
long-range one-pion-exchange one. A possible way to improve the situation would be to 
calculate higher order corrections employing the modified counting rules or to employ a 
relativistic formalism which we will discuss later. Because performing higher order loop 
corrections is a formidable task, an easy way to circumvent the problem might be to 
employ a pionless theory in which virtual pions exchanged between the two nucleons are 
integrated out; in this pionless theory, the one-pion exchange, two-pion exchange and 
contact terms in HB^PT are subsumed in a contact term. 

In this work we employ a pionless effective field theory (EFT) with a di-baryon [27, 28]^ 
and external pions [31] to calculate the total cross section of the pp — > ppn'^ process. 
A main motivation of this work stems from the observation that some of the HB^PT 
calculations failed in reproducing the experimental data. Meanwhile, it is known that 
the energy dependence of the experimental data can be well reproduced in terms of the 
final state interaction and the phase space [2]. The pionless theory would be a "minimal" 
formalism to take account of these two features. Furthermore, after taking these two 
features into account, the difference between the theory and experiment appears in the 
overall factor and the experimental data can be easily reproduced by fitting an unknown 
constant that appears in a contact vertex. 

This paper is organized as follows. In Sec. 2 the pionless effective Lagrangian with a 
di-baryon and external pions is introduced. In Sec. 3 we fix low-energy constants (LECs) 
for the final ^5*0 and initial ^Pq NN states. In Sec. 4 the amplitudes of pp — > pp7r° for a 
leading one-body and subleading contact vertex are obtained including the strong initial 
state interaction and the strong and Coulomb final state interactions. We estimate in 
Sec. 5 the value of an LEG in the external neutral pion production contact vertex from 
HBxPT, and in Sec. 6 we show our numerical results of the total cross section and compare 
with the experimental data. Finally, in Sec. 7, the discussion and conclusions of this work 
are given. 

2. Pionless effective Lagrangian with a di-baryon and external pions 

^Some other issues, a role of three-particle singularities for the ttNN system (using a toy model) [25] 
and the use of the effective low-momentum NN potential Vioiv-k for the process [26]. were also studied. 

^We have studied np dj cross section at BBN energies [29] and neutron-neutron fusion process [30] 
employing this formalism. 
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An effective Lagrangian witliout virtual pions and witli a di-baryon and external pions 
for describing the pp — > ppir^ reaction may read 



£, — Cm + Cs-\- Cnstt + ^ 



(1) 



where >Cjv is the one-nucleon Lagrangian interacting with the external pions, Cs is that 
for the ^5*0 channel di-baryon field interacting with two-nucleon, Cnsw is the contact in- 
teraction Lagrangian for external pion-dibaryon- two-nucleon, and C^^^ is the two-nucleon 
Lagrangian for the ^Pq channel. 

The one-nucleon Lagrangian £n in heavy-baryon formalism reads 

Cn ^N^ {iv-D + 2igAS ■ A + \{v ■ Df - + 2gA{v ■ A, S ■ D} + ■■■]] N, (2) 

where is a velocity vector with a condition v'^ = 1 where = (1,0), and the spin 
operator is 5^ = (0,a/2). D^ = d^ + with = and A^ = 

The nonlinearly realized external pions in the ^ field are given by ^ = exp[ir • 7r/(2/7r)]. 
gA is the axial vector coupling. 

The effective Lagrangian for the two-nucleon part may read [27, 28, 32, 33] 



1 



ivD + [{v ■ Df - D^] + 5s 



Am 



N 



4(7V^Pi'^o)^)+h.c.l,(3) 



d^3 



■'NsTT — 



N^a2a ■ i{D - D)T2nN 



NN 



VSmTvro 

c%Ai-An^oZn)\n-oZn) + 



[iv ■ Ac) -|- h.c. 



with 



O 



ij,a 



i(r). p(P) _ p^"^) r,.) 



(P) 



Vs 



(4) 
(5) 

(6) 



where Sa is the (spin singlet) ^5*0 channel di-baryon field with the isospin index a, ag is 
the sign factor = ±1. is the mass difference between the di-baryon mass and 
two-nucleon mass, i.e., = 2mN -\- Sg- Us is the coupling constant of the di-baryon and 
two-nucleon interaction. P/*^^ are the projection operators for S — ^Sq and ^Pq channels; 

Pi^-o) ^ 1^ , p(3Po) ^ 1 . , / ^ E Tr (p(-)tp(-')) ^ U^s' 

V o V O pol.avg 

where ai (r^) is the spin (isospin) operator, df^ is an unknown LEG of the (external) 
pion-(spin singlet) dibaryon-nucleon-nucleon {nsNN) interaction. Tq is the effective range 
in the ^5*0 (pp) channel, and A'^ = Y^a- ^2 is LEG for the P-wave A^A^ scattering 
in the ^Pq channel. 



3. Fixing LECs of the initial and final NN interactions 
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Figure 1: Diagrams for the dressed di-baryon propagator including the Coulomb inter- 
action. A double-line with a filled circle denotes the renormalized dressed di-baryon 
propagator. Double-lines without a filled circle and single-line-curves denote the bare di- 
baryon propagators and nucleon propagators, respectively. The two-nucleon propagator 
with a shaded blob denotes the Green's function including the Coulomb potential. A 
(spin-singlet) dibaryon-nucleon-nucleon (sNN) vertex is proportional to the LEC i/s- 



In this section, we calculate the S- and P-wave NN scattering amplitudes to fix the 
LECs in the two-nucleon part. In Fig. 1, diagrams for the dressed ^5*0 channel di-baryon 
propagator arc shown where the two-nucleon bubble diagrams including the Coulomb 
interaction are summed up to the infinite order. The inverse of the propagator in the 
center of mass (CM) frame is given by 

iD;\p) = ias{E + 5s) - iyUoip) : (8) 

with 



where G^^^ is the outgoing two-nucleon Green's function including the Coulomb potential, 

G^c\E)^ ^. , (10) 

where E is the total CM energy, E = p^/mjv + ■ ■ ■, Hq is the free Hamiltonian for two- 
proton, Hq = p'^/ruN, and Vc is the repulsive Coulomb force Vc = a/r: a is the fine 
structure constant. Employing the dimensional regularization in d — A — 2e space-time 
dimensions, we obtain [34, 35] 



Mp) = 



am 



Stt 



1 , / TT/i^ 

37 + 2 + In ' ^ 



-^Kv)-c,—{ip), (11) 



where // is the scale of the dimensional regularization, 7 = 0.5772 • • •, and 

00 -j^ 

h(r]) ^ Rei/:{ir)) -Inrj, Reip{r)) = rfY^ , ^ 2\ ~^^ 

v=l + ^7 J 

S-e2-._l' 2p ■ ^^^^ 
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Figure 2: Diagram for the S'-wave pp scattering amplitude with the Coulomb interaction. 
See the caption of Fig. 1 for details. 



Thus the inverse of the renormalized dressed di-baryon propagator is obtained as 



iD:\p) 



rriNy; 



+ 



^ „2 



p + amNh{r]) + ip 



(13) 



where 5f is the renormalized mass difference between the di-baryon and two nucleons, 



R 



,am 



N 



Stt 



1 



- 37 + 2 + In 



TT/i 
2^2 



N , 



(14) 



We fix it by using the scattering length ac below. 

In Fig. 2, a diagram of the 5*- wave pp scattering amplitude with the Coulomb interaction 
is shown and thus we have the ^-wave scattering amplitude as 



iAs = {-iysi^o)iDs{p){-iysi^Q) 

. An 



— i- 



fi2 2iao 



rriN 



rriNVs 



— 1 — 9" 



amNh{r]) — ipC^ 



(15) 



with 



(16) 



where {k\ilj^^) arc the Coulomb wave functions obtained by solving the Schrodinger equa- 
tion {H — E)\'iIj^^^) = with H = Hq + Vc and represented in the \k) space for the two 
protons. (Jo is the yS-wave Coulomb phase shift (Tq = arg r{l + irj). The 5"- wave amplitude 
As is given in terms of the effective range parameters as 



47r 



iT^N--^ + \rop^ H aniNhivi) -ipC^' 



(17) 



where ac is the scattering length, ac = —7.8063 ±0.0026 fm, tq is the effective range, tq = 
2.794 ± 0.014 fm, and the ellipsis represents the higher order effective range corrections. 
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Figure 3: Diagrams for the P-wave NN scattering. Two-nucleon bubble diagrams are 
summed up to the infinite order. Lines and curves denote nucleon propagators. A four- 
nucleon contact vertex with a filled box is proportional to the LEG Cg. 



Now it is easy to match the LECs with the effective range parameters. Thus we have 
(Ts = —1 and 

, 2 /27r mNTo 1 , . 

TUN \ To 2 — - 2^0^^ + amNh{ri) + ip 

We note that the sign of the LEG ys cannot be determined by the effective range param- 
eters. 

In Fig. 3, diagrams for the P-wave NN scattering are shown. Because the momenta 
of the two protons are quite large for the pion production reaction, we cannot treat the 
P-wave vertex correction in a perturbative way. Thus the two-proton bubble diagrams 
are summed up to the infinite order without including the Goulomb interaction ^ and 
the LEG 6*2 is renormalized by a phase shift at the threshold energy. The scattering 
amplitude for the ^Pq channel is obtained as 

- ip' 

where we have used the nucleon propagator, iSN{k) — i/[ko — k'^/{2mN) + ^e], where 
is the residual nucleon momentum k'^ — — rriNV^. is the nucleon momentum 
P^ = m^, and we have employed the dimensional regularization for the loop calculation. 
The LEG C2 is fixed by the phase shift of ^Pq channel at pion production threshold, 
^pipth) — —7.5° at pth — \/'rfi-KTajq . Thus we have 

47r 

p,\cot5p(pi^). (20) 







m]siC\ 

4. Amplitudes for pp ppn^ near threshold 

In Figs. 4 and 5, we show diagrams for pp — > ppir^ near threshold. In diagram (a) in 
Fig. 4 and (c) in Fig. 5, the pion is emitted from the one-body ttNN vertex. Because 
the S'-wave pion production is considered, the ttNN vertex (the vertex with "X" in the 
figures) is obtained from the l/rriN Lagrangian. In the diagram (b) in Fig. 4 and (d) in 



""■An effective Coulomb interaction is denoted by 77 = aM/{2p). For the high momenta p > aM/2 ~ 
3.42 MeV, 77 < 1 and the Coulomb interaction can be treated in a perturbative way. 
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P P 

(a) (b) 

Figure 4: Diagrams for pp — > ppn'^ near threshold with the strong and Coulomb final state 
interactions and without the initial state interaction. 5'-wave neutral pion is emitted from 
nNN vertex with "X" at the I/uin order in (a), whereas the pion is emitted from a pion- 
dibaryon-nucleon-nucleon contact vertex in (b) which is proportional to LEG 



(c) (d) 

Figure 5: Diagrams for pp — > pp7r° with the strong initial and the strong and Coulomb 
final state interactions. See the caption of Fig. 4 for more details. 




Fig. 5, the pion is emitted from the pion-(spin singlet) dibaryon-nucleon-nucleon (ttsNN) 
contact vertex which is proportional to the unknown LEC d^^ . The strong and Coulomb 
final state interactions are considered in all of the diagrams (a), (b), (c) and (d) in Figs. 4 
and 5, whereas the strong initial state interaction is considered ^ in the diagrams (c) and 
(d) in Fig. 5. 

The one-body amplitude from the (a) and (c) diagrams and the two-body (contact) 
amplitude from the (b) and (d) diagrams are obtained as 



^^n!^ ^ - \rQV''^ + amNh{ri') + ip' C^, i _ ^^^ip 

■A =4 f^SL C^'e^^^uj^p 1__ 

^-^ib+d) y _L _ Irop'^ + arriNhirj') + ip' C^, i _ rn^i 



(22) 



where we have used the nucleon propagator, iS^ik) = i/[kQ — k'^/{2mN) + ie] in the 
calculation of the (a) , (c) , (d) diagrams. 2p and 2pl are the relative three momenta between 
incoming and outgoing two protons, respectively; p — \p\ and = Ip*]. r]' — aniN/i^p') 



^See the footnote 4. 
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Figure 6: Diagram for a one-pion exchange contribution to the pp — > pjm^ process for 
estimation the LEG in the contact vertex. A double-hne, single-hne, and dashed-hne 
denote a di-baryon, nucleon, and pion, respectively. A irNN vertex with a dot is obtained 
from C^^lf and a tttiNN vertex with "X" is from C^^jf and Cflf. 



and Uq is the energy of outgoing pion, Uq — y^g^ + m^: q is the outgoing pion momentum. 
In the loop calculations for the diagrams (c) and (d) in Fig. 5, we have employed the 
dimensional regularization, and we have neglected terms involving q and p', which are 
considered to be small as compared to p ~ yfWi~m^. We note that there remain no 
unknown parameters in the amplitudes except for the LEG d^^'> in the two-body (contact) 
amplitude in Eq. (22). 

5. Estimate of the LEG from HBxPT 

In this section we estimate an order of magnitude of the LEG d^^^ from HBxPT. We here 
consider a one-pion-exchange (OPE) diagram shown in Fig. 6. This diagram is the lowest 
order OPE contribution in the standard Weinberg counting rules. We include a higher 
order (relativistic) correction to the ttttNN vertex which is found to be important [20] and 
is, in the modified counting rules, of the same order as the lowest order diagram. Because 
other diagrams, e.g, two-pion exchange diagrams would give comparable contributions to 
that of the LEG 4^) 

(see Refs. [17, 18]), we expect that the estimate of the diagram in 
Fig. 6 would be reliable only for the order of magnitude estimation for the LEG 
The effective chiral Lagrangian relevant to this purpose reads 



(23) 



where is the xPT Lagrangian for the pions, and Ct^n is the HBxPT Lagrangian for 
the pions and nucleon. These Lagrangians are expanded as 

£. = 4^) + ---, £.^, = £a + 4'ii + 4'i + ---, (24) 

where is the standard chiral Lagrangian for pions at LO. The expression of C^^lf and 
some part of £j|.jv given in Eq. (2). To calculate the isoscalar tttiNN interaction in 
the diagram in Fig. 6, we consider the interaction Lagrangian [36] 



,(2) 



,(3) 



= Art 



ciTrix+) + 



9l 



4:171% 



2m]si 
iv ■ AA • D + h.c. 



- 4c2 j {v ■ A)2 - 4C3A • A 

2C2 



TUN 



[iTr(v ■ AA^)D^ + h.c] 



N + 



(25) 
,(26) 
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(3) 

where we have included rclativistic corrections in the higher order Lagrangian C)^^ as 
mentioned above. The values of the LECs Ci , C2 and C3 are fixed in the tree-level calcu- 
lations [37] ^ as 

ci = -0.64, C2 = 1.79, C3 = -3.90 [GeV-^] . (27) 
Thus a contribution to the ttsNN vertex function from the diagram in Fig. 6 is obtained 



as 



.-^ _ .mNysgA\p\eabc 

TTsNN — —1'- 



-4ci + 2c2 - — ^ + cs]ml dx- 
IbiriN 



(28) 
with 



167r/3 



F ^ x{l- x)p^ + xmi, (29) 

where a and b (in Cabc) are the isospin indices for the initial and final two-nucleon state, 
respectively, and c is that for the outgoing pion. In the calculation of the diagram in 
Fig. 6, wc have used the nuclcon propagator, iS]\f{k) = i/[kQ — k^/{2mi\f) + ie], and the 
"potential" pion propagator, iA^(A;) = —i/{k'^ -|-m^ ^^^) because the typical momentum 
and energy transfer between two protons are P ~ m^mAr and ko ~ m7r/2, respectively. 
Furthermore because ~ rriT^mN » iri^, we take an approximation F ~ a:(l — x)p^ and 
thus have 

Therefore the value of the LEG d'^^ from the loop diagram in Fig. 6 is obtained as 



nsNN - Tir^f^^ ( -4ci + 2C2 - + C3 eabc ■ (30) 



where the different signs for d^^ have been obtained because of the LEG ys in Eq. (18). 

6. Numerical results 

Total cross section of pp — > pp7r° near threshold is calculated using the formula 

1 yg-""^ da da 1 m^YgV >r- , .,2 /qo^ 



spin 



with 



P 



— mE 



T ' 



(33) 



2mjv 



^Values of the LECs ci, C2, and C3 fixed in the one- loop calculations are quite different from those in 
the tree-level ones. See, e.g., Refs. [11, 38] for details. 
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where q is the outgoing pion momentum, q= 2p and 2p' are relative momenta of the 
initial and final two protons, respectively. T is the initial total energy T ~ /rriN and 
viab — 2p/m;v- We have expanded the proton energies in the phase factor in terms of 
l/rriN and kept up to the l/rriiy order. A is the amplitude A = A(a+c) + -^(b+d) where 
-^(a+c) and A(^b+d) arc obtained in Eqs. (21) and (22), respectively. Note that the factor 
1/2 in Eq. (32) is the symmetry factor for the final two protons. 

In Fig. 7 we plot our results of the total cross section as a function of rj^^ — g"*"^/m^. 
The solid curve and long-dashed curve have been obtained by using d^^ — ±0.140 fm^/^ 
fixed from the one-pion exchange diagram in Fig. 6 in the previous section. The LEG 4^) 
is also fixed by using the experimental data as 

di^)med ^ _Q j2, +o.55fm^/^ (34) 

where we have two values of d!^^ with different signs. The short-dashed curve is obtained 
by using o^(^2)/i«ed _ _q 12 fm^/^. The dotted hue corresponds to the case where only the 
contribution from the one-body amplitude A^a+c) is considered. The experimental data 
are also included in the figure. 

Wc find that the experimental data are reproduced reasonably well with the value of 
d^^ = —0.14 fm^/^. By contrast, we obtain almost vanishing total cross sections with the 
value d^^^ = +0.140 fm^/^ because the two-body amplitude with d^^^ — +0.140 fm^/^ is 
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almost canceled with the amplitude from the one-body contribution. On the other hand, 
for the whole energy range the experimental near-threshold cross section data are well 
reproduced with the use of the fitted parameter c^(^2)/i«erf _ _q 12 fm^/^. 

Wc also find that approximately a half of the experimental data comes from the one- 
body (lA) amplitude in the pionlcss theory. Whereas, in the previous DWBA calculations, 
about 1/5 of the experimental data comes from the one-body matrix element. The dif- 
ferent one-body contributions in the pionful and pionless theory stem from the different 
short range contributions in the pionless and pionful theory. The importance of the short 
range contributions in the one-body part of the DWBA calculations has been pointed out 
(see, e.g., in Fig. 10 in Ref. [13]). In the pionless theory, we conjecture that some short 
range part of the one-body matrix element in the DWBA calculations is integrated out 
and held in the contact LEG . 

7. Discussion and conclusions 

In this work we calculated the total cross section for pp — > pp7r° near threshold in 
pionlcss EFT with the di-baryon and external pion fields. The leading one-body ampli- 
tude and subleading contact amplitude were obtained including the strong initial state 
interaction and the strong and Coulomb final-state interactions. After we fix the LECs 
for the A'^A'^ scatterings, there remains only one unknown constant, d!^\ in the amplitude. 
We estimated it from the one-pion exchange diagram in the pionful theory. Although 
this method does not allow us to fix the sign of we have found that one of the two 
choices for d^^ leads to the cross sections that agree with the experimental data reasonably 
well. On the other hand, the whole range of the experimental data near threshold can be 
reproduced by adjusting the only unknown LEG in the theory, d^\ As discussed in In- 
troduction, this is an expected result because the energy dependence of the experimental 
total cross section is known to be well described by the phase factor and the final-state 
interaction [2], which have been taken into account in this work, and the overall strength 
of the cross section can be adjusted by the value of d'^\ This feature would be the same 
in the NNLO HB^PT calculations because an unknown constant appears in the contact 
ttNNNN vertex and can be adjusted so as to reproduce the experimental data (though 
there are many other corrections coming out of the pion loop diagrams). 

The main problem in the NNLO HB^PT calculations that the contributions of formally 
different chiral orders become of comparable magnitude has not been fully clarified in this 
work. Though the main difficulties in the HBxPT calculations may stem from the sup- 
pression mechanisms in the production operator, another difficulty would stem from the 
heavy-baryon formalism involving the typical large momentum transfer. The equations 
of motion of the heavy-field, t; ■ p ~ fp' / {2m]\j), makes a connection between the terms in 
the different orders, one derivative term v -p and two-derivative term p^/[2mN), and thus 
the order counting rules become not transparent.^ This is the reason why we needed to 
take account of the relativistic corrections: e.g., the term proportional to v ■ A A ■ D in 
in Eq. (26) is a relativistic correction to the term proportional to (v ■ A)"^ in in 

''When there is one heavy-particle, we can avoid the problem by choosing the velocity vector so as 
to V ■ p = 0. But if there are two heavy-particles, we cannot avoid the problem. 



12 



Eq. (25). As mentioned in Introduction, one way to solve the problem will be to employ 
the modified counting rules (see Refs. [23, 24] as well) and calculate the cross section 
collecting all of the pieces up to NNLO. However, because a lot of higher order terms 
will be involved in the NNLO HB^PT calculation with the modified counting rules, an- 
other possible way to calculate the production operator may be to employ a manifestly 
Lorentz invariant baryon chiral perturbation theory [39] with an additional subtraction 
scheme, such as the infrared renormalization scheme [40] or the extended on-mass shell 
scheme [11]. That would be worth studying in order to clarify the issues pertaining to 
the HBxPT calculation of the near-threshold pp — > ppir^ reaction up to one-loop order. 
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